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SUMMARY 

Une approche pour calculer l'ecoulement uniforme et les characteri stiques 
dynamique des fermetures etanchs et des paliers est discuter. La velocitee 
moyenne locale est tres affectais par les effets d'entree et de sortie et 
1 'injection liquide, ce qui cause des ecoulements secondaires. Pour les 
maquettes en trois dimension considerez, les resultats integrales moyens etait 
en accord avec les resultats d'analyse. Les mesurements de pression en 
direction uniques n'etait pas suffisant pour donnes les variations 
d ' ecoulement. Mai, pour les correlations de fuites d'etanches et de paliers, 
les principes d'etats de correspondence etait utiles. Les trois phenomenes 
trouvait pendant les epreuves de trois configurations de fermeture etanche 
excentrique sans| rotation pour le "Space Shuttle Main Engine" (SSME) programme 
sont aussi discutez. L'injection des fluides, callages de fermeture etanche, 
et profile de pression sont examiner en relation au zones d'ecoulement 
secondaires de separation et de rigidite. 

1. INTRODUCTION 

Although the primary function of a seal is to control leakage, a secondary 
but equally Important purpose is to provide dynamic stability. A seal or 
bearing properly designed by using computational fluid dynamics and knowledge- 
based methodologies can enhance turbomachine stability. 

Much of the current research on seal- or bearing-driven Instabilities in 
turbomachines is contained in NASA CP's 2133, 2238, 2250, 2443, and 2409 
(refs. 1 and 2). Topics such as seal dynamic coefficients (as measured by 
Childs, Benckert, and lwatsubo), analytical methods (Fleming, Nelson, Black, 
Brown, Muszynska, and Bently), numerical techniques (CHAM, Rhode, and 
Nordmann), and practical applications (gas injection and swirl brakes) are but 
a sampling of the research efforts of the past decade. Leakage rates for 
cylindrical and tapered seals can be predicted reasonably well, but direct 
stiffness predictions are not very accurate as revealed by a comparison with 
noncryogenic data (ref. 3). Tests with nonrotating convergent tapered seals 
revealed a rapid dropoff in direct stiffness for Inlet/exit area ratios greater 
than 2 (ref. 4). Parallel or series hydrostatic bearings and magnetic bearings 



appear to have potential for high DN and dynamic applications, but stability 
and longevity remain as problems (refs. 1, 2, and 5). 

Herein we outline a relativistic method to calculate flow and dynamic 
characteristics for seals or bearings and compare this solution method with 
selected data of Childs (ref. 3). We also discuss leakage rates and relations 
between pressure profile data for three eccentric nonrotating seal tests, 
associated configurations, and corresponding states. 

2. NUMERICAL APPROACH 

Consider a shaft at a static eccentricity a and a dynamic eccentricity 
d rotating clockwise within a static housing (fig. 1). Such a configuration 
can represent a seal or a bearing. The axial pressure drop for the bearing Is 
quite small, but for the seal It Is large and engenders significant Coriolis 
and centrifugal forces at the Inlet or convergent/divergent zones. 


The problem here Is that we are dealing with nonconventlonal flow fields 
In very narrow passages with the passage geometry changing periodically with 
time. But In order to simulate seal or bearing flows, the coordinate system 
must be transformed Into the relativistic frame such that the geometry is 
invarlent in time (refs. 6 and 7). 

To an observer in the rotating frame the relation between absolute, rela- 
tive, and grid velocity becomes 
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where eq. (5) is for constant w. The resulting equation of motion becomes 
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Expanding this form of the Navier-Stokes equations gives rise to two source 
terms to be added to the existing cylindrical polar numerical code. These 
source terms for the U and V momentum equations become 

S(U) = -2pwVVol ce ii CO 

S (V) = (pco 2 r + 2pwU)Vol ce il (8) 

The turbulent flow in a seal is simulated by way of the Prandtl mixing- 
length model. The mixing length a. m is characterized as 

8 . 

-r- = K * K = 0.435 

o o 
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where <S is the normal distance measured between the stator and the rotor 
surface and y is the minimum local distance measured from the stator and the 
rotor surfaces. 


The resulting Navier Stokes equations describe the time-averaged 
distributions of velocity, pressure, etc. Note that the str ess t ensor 
described earlier now also contains Reynolds stress terms -pu^. The 

turbulence model currently used in our seal flow simulation makes use of the 
eddy- vi scosity concept to compute the Reynolds stresses from 


-pu'uj = y t 




(10) 


where is the Kronecker operator, y-t = pS-fn'P the ec| dy viscosity with 

8, m the mixing length, and ip is the rate of strain. 


If w e de fine the effective viscosity as yeff = y| + vt» the stress tensor 
- pu^uj) can be simply replaced by V(-r e ff). The shaft reaction forces 
are 


Radial : 


Tangential : 
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Further, In the numerical approach the shaft rotates clockwise, and the tan- 
gential force F e would be defined by -F e . Stiffness is then related to the 
ratio Fr/clearance, and damping is related to the slope of the Fe/clearance- 
versus-rotational-speed locus. The application of these equations to seal data 
is discussed later. For purposes of the discussion u = Q. 

3. DISCUSSION 

Figures 2 and 3 present the distribution of radial and tangential forces 
as a function of rotational speed for the experimental data of Childs (ref. 3) 
and that calculated numerically. Although there were three steps within the 
length of the experimental seal and the modeled seal was considered as straight 
cylindrical, without inlet preswirl, experiment and calculations agreed rea- 
sonably well. The calculated radial forces Increased with speed but did not 
exhibit the positive values of the experimental data at elevated rotor speed. 
The omission of the internal steps or preswirl may account for this discrep- 
ancy. The tangential forces -F e Increased with rotational speed and agreed 
reasonably well with the data of Childs (ref. 3). 

The discussion now centers on three phenomena associated with three dif- 
ferent fully eccentric nonrotating seal configurations. 

First, helium was Injected into a nonrotating straight cylindrical seal 
at the -90°, 0°, and 90° positions at approximately 0.8 of the seal length 
(fig. 4(a)). The overall magnitude of the pressure over the surface at the 0° 
position Increased, but the magnitude over the surface at 180° decreased 
(fig. 4(b)). Such profiles are similar to those for an eccentric rotor; thus 
the Injected fluid essentially functioned as an eccentric rotor. We think that 
the velocity patterns moved in a circumferential manner and that for fluids the 
Injection engenders two-phase flow. Such patterns would indicate that injec- 
tion direction is critical. Injection opposite to the direction of rotation 
would enhance stability; injection in the direction of rotation would be 
destabilizing. This dynamic stabilizing technique could be applied to bearings 
or seals (refs. 1 and 2). 

The second phenomenon was found in the second step of a nonrotating three- 
step seal (fig. 5(a)) where the third step was choked in the 180° position. 
Within the first step the pressure profile at the 0° position was larger than 
that for the 180° position. However, within the second step the profiles 
crossed over (fig. 5(b)). We think that the flow at the 180° position was 
partially blocked by the two-phase choke at the entry to the third step and was 
diverted toward the 0° position. The flow was then rediverted within the third 
step toward the 180° position, and the pressure profile mimicked that of a 
diffuser (i.e., there was a monotone Incr e a se of pressure with axial position). 

Such fluid jetting was also found to exist in single and sequential ori- 
fice and Borda Inlets. It was assumed that liquidlike behavior produced fluid 
jetting but that data exhibiting gaslike behavior were produced by the strong 
separation effects (refs. 8 and 9). Fluid velocities for liquids with high 
subcooling fluid velocities can exceed the local sonic velocity and remain 
"supersonic 11 to 105 L/D. With proper care a seal could be designed to provide 
for such separations and thus enhance both stability and leakage control. 
However, the effect of shaft rotation on jetting zones is yet to be determined. 
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The third phenomenon was found within the steps of a three-step labyrinth 
seal of 12, 11, and 10 teeth, respectively, (fig. 6(a)). The carryover 
increased with axial position; so although the pressure difference between 0° 
and 180° was positive at the inlet, it became negative at the exit for each 
step (fig. 6(b)). Such crossovers in the pressure profiles provided positive 
and negative direct stiffness within each step with little or no resultant 
stiffness. 

It is evident from previous discussion that fluid conditions at the inlet 
and exit can create zones of secondary flows that either enhance or diminish 
the average fluid velocity. Further, for the eccentric nonrotating tests dis- 
cussed herein the pressure measurements were taken at 0° and 180° along the 
axial length, but these positions may not represent maximimum-to-mlnlmum pres- 
sure variations. For example, when eccentric seals or bearings rotate at high 
speeds, the maximum and minimum pressures may be quite close geometrically. 

This points out the need for careful instrumentation guided by three- 
dimensional computations, data, and intuition. 

4. KEY TO STABILITY 

The key to stability is the average velocity profile. As the local aver- 
age velocity within the bearing or seal decreases, the cross-coupled terms that 
contribute to Instability are decreased (ref. 10). However, the Increase in 
damping that is associated with decreases in the cross-coupled terms usually 
effects a decrease in direct stiffness. Thus one must consider system dynamics 
and not just those of the seal or bearing alone. For a simulated bearing 
Prof. M.O . Braun of the University of Akron, Akron, Ohio, in a private commun- 
ication has demonstrated that zones of secondary flow exist in the inlet and 
exit regions of a bearing such that the average fluid velocity is a minimum 
near the minimum gap. 

From these results it is apparent that, although some average set of 
dynamic coefficients may be found, the three-dimensional effects are signifi- 
cant and local variations in these coefficients must be properly factored into 
an analysis. 


5. CORRESPONDING STATES 

This final section could also be entitled "similitudes that are useful to 
design and guides for experimentation." In most high-performance turbomachine 
applications the pressure gradients are sufficiently large to achieve sonic 
velocities (refs. 1 and 2). The principles of corresponding states have been 
applied to single- and multicomponent fluids and to single- and two-phase 
choked flows through venturis and through single and sequential Borda and 
orifice configurations (refs. 1, 9, 11, and 12). The principles have also been 
applied to regions of annular flows as found for example in bearings and seals 
(refs. 4 and 13 to 16). Reduced mass flux data for three space shuttle main 
engine (SSME) seal configurations for the cryogens nitrogen and hydrogen 
(refs. 13 to 15) tend to verify extending the principles of corresponding 
states to other fluids such as oxygen. However, for the extensions to Include 
gaseous helium, the reduced helium pressure must be further normalized by a 
factor of TO (fig. 7). At this time the reason for the factor of 10 is 
undetermined. 
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6. SUMMARY OF RESULTS 


Herein we have discussed an approach to calculating flow and dynamic 
characteristics for seals or bearings. For the restricted model considered, 
the results are In reasonable agreement with a run selected from the data of 
Childs (ref. 3). 

Me discussed leakage rates and relations between pressure profiles and 
direct stiffness for three seal configurations of the space shuttle main engine 
class and related these profiles to postulated zones of secondary flow or sep- 
aration. For these tests the seals were eccentric but the shaft was not 
rotating. Fluid Injection was found to augment direct stiffness. Flow choking 
can occur within a seal. Within stepped configurations choking diverted the 
flow and effected a crossover In the pressure profiles within the upstream 
step, thereby decreasing direct stiffness (the same comment applies to the 
straight cylindrical seal). For a labyrinth seal carryover diminished from the 
entrance and effected a crossover In the pressure profiles within each of the 
three steps. As such, the direct stiffness changed from positive to negative. 

Pressure measurements taken at 0° and 180° along the axial length of a 
seal or bearing are Insufficient to define maxlmlmum- to- minimum pressure vari- 
ations and zones of secondary flows. This points out the need for careful 
Instrumentation guided by three-dimensional computations, data, and Intuition. 

The local average velocity within a bearing or a seal emerged as the key 
to stability and was strongly Influenced by Inlet and exit effects and fluid 
Injection, which In turn drove zones of secondary flow. Thus although some 
average set of dynamic coefficients may be found, the three-dimensional effects 
were significant and local variations In these coefficients must be properly 
Integrated Into an analysis. 

The principles of corresponding states have been applied to seal and 
bearing flows and were found to be useful. 
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Figure 1. - Schematic of simulated stepped-seal configuration 
in a constant-Z plane. 


AXIAL 

REYNOLDS 

NUMBER, 

lxlO 3 r A 



ROTOR SPEED, u, rad/sec 

Figure 2. - Radial force as function of rotor speed for two axial Reynolds numbers. 
Stepped-seal data of Childs (ref. 3). (Herein u = Q. ) 
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Figure 3. - Tangential force as function of rotor speed for two axial Reynolds 
numbers. Stepped-seal data of Childs (ref. 3). Whirl frequency ratio, 

k YU /Cg = as u increases. (Herein u = Q ) 
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(b) Pressure profiles with and without fluid injection, 
Figure 4. - Fluid injection into straight cylindrical seal. 
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(a) 33-Tooth labyrinth shaft seal configuration. (Dimensions are in centimeters.) 
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(b) Pressure profiles for gaseous hydrogen with crossover; 2/3 fully eccentric position 
Figure 6. - Three-step labyrinth seal. 
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Figure 7. - Reduced mass flux as a function of reduced pressure differential for gases nitrogen, hydrogen, 
and helium flowing through a straight cylindrical seal (see fig. 4). 
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